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Abstract 


A model  for  prediction  of  burning  rate-pres- 
sure-crossflow velocity  relationships  for  non-met- 
alized  composite  propellants  containing  unimodal 
oxidizer,  given  only  composition  and  oxidizer  par- 
ticle size,  has  been  developed.  This  model  embodies 
many  of  the  concepts  used  in  the  Becks tead-Derr- 
Price  model,  but  contains  major  modifications,  in- 
cluding a postulated  columnar  diffusion  flame  bend- 
ing mechanism  for  erosive  burning.  The  major  part 
of  this  paper  is  devoted  to  description  and  discus- 
sion of  these  modifications  and  to  model  development. 
Preliminary  predictions  of  burning  rate  at  various 
pressures  and  crossflow  velocities  have  been  made 
for  a series  of  three  73/27  ajasonium  perchlorate 
(AP) /hydroxy-terminated-polybutadiene  (HTPB)  formu- 
lations, with  oxidizer  particle  diameters  of  5,  20, 
and  200  microns , and  compared  with  data  for  these 
formulations.  With  optimization  of  three  ♦free 
constants'*  appearing  in  the  model,  it  is  found  to 
give  excellent  agreement  with  no-crossflow  burning 
rate  data  over  the  entire  range  of  pressures  and 
particle  sizes  studied.  In  all  cases,  however,  the 
predicted  sensitivity  of  burning  rate  to  crossflow 
velocity  is  somewhat  less  than  observed  experiment- 
ally. 


Introduction  and  Background 


As  part  of  a program  in  which  the  sensitivity 
of  composite  propellant  burning  rate  to  crossflows 
is  being  experimentally  and  analytically  studied, 
this  author  has  developed  a "first  generation" 
model*  for  the  prediction  of  burning  rate-pressure- 
crossflow  relationships  for  composite  propellants 
which  requires  as  input  burning  rate  versus  pressure 
data  in  the  absence  of  crossflow.  A more  fundamen- 
tal model  (with  explicit  calculation  of  the  dis- 
tances of  various  heat  release  zones  from  the  pro- 
pellant surface  rather  than  inference  of  these  dis- 
tances from  zero  crossflow  data)  of  the  propellant 
combustion  process  which  permits  prediction  of  burn- 
ing rate  versus  pressure  at  various  crossflow  vel- 
ocities (including  zero  crossflow) , given  only  pro- 
pellant composition  and  ingredient  size  data,  is 
highly  desirable.  During  the  past  year,  this  author 
has  been  working  on  development  of  such  a model  for 
the  case  of  propellants  containing  unimodal  ammon- 
ium perchlorate  oxidizer  and  no  metal  additives 
(with  plans  to  extend  this  model  later  to  treat  mul- 
timodal oxidizer  and  metallized  propellant  cases). 
Initially,  it  was  planned  to  simply  modify  the 
Beckstead-Derr-Price  (BDP)  model2*3  (the  most  com- 
monly accepted  composite  propellant  combustion  model, 
in  various  forms,  in  this  country  today)  for  pre- 
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diction  of  burning  rates  of  composite  propellants 
as  a function  of  pressure  (in  the  absence  of  cross- 
flow)  to  allow  for  bending  of  the  diffusion  flame(s) 
considered  in  that  model.  However,  upon  careful 
review  of  the  BDP  model,  this  author  found  sufficient 
problems  and  areas  of  disagreement  with  that  model 
that  it  was  decided  to  develop  an  entirely  new  com- 
posite propellant  combustion  computer  code  (embody- 
ing many  of  the  BDP  concepts,  while  modifying  or 
replacing  others)  with  the  flame-bending  mechanism 
described  in  Reference  1 embedded  in  the  mathemat- 
ical analysis.  Major  modifications  to  the  BDP  model 
included  are: 

1)  Variation  in  local  oxidizer/binder  surface 
area  ratio  as  the  propellant  regresses  past 
an  oxidizer  particle  is  considered.  (In  the 
BDP,  an  average  ratio  is  used  - this  assumes 
that  several  very  nonlinear  processes  can  be 
linearly  averaged.) 

2)  The  kinetics  of  subsurface/surface  exothermic 
reactions  are  considered,  with  use  of  rate 
expressions  based  upon  the  work  of  Waesche 
and  Wenograd  . (In  the  BDP  model,  subsur- 
face/surface heat  release  is  included  with 
the  endothermic  ingredient  vaporization  heats, 
with  the  resultant  implicit  assumption  that 
the  amount  of  heat  release  in  these  reactions 
per  unit  mass  of  propellant  is  Independent 

of  such  parameters  as  burning  rate.) 

3)  A correction  of  an  inconsistency  in  defini- 
tion of  areas  in  the  BDP  model  is  made. 

4)  The  calculation  of  the  dimensionless  stoich- 
iometric group  needed  for  calculation  of  the 
flame  height  via  the  Burke-Schumann5  analysis 
is  modified.  (The  group  used  in  the  BDP 
model  is  inconsistent  with  that  defined  in 
the  original  work  of  Burke-Schumann.) 

5)  A two-flame  (fuel-gas/oxidizer-gas  columnar 
diffusion  flame  and  ammonium  perchlorate  mono- 
propellant flame) , rather  than  a three-flame 
model,  is  used.  (With  correction  of  the  cal- 
culation of  the  stoichiometry  dimensionless 
group  for  the  Burke-Schumann  analysis,  it  no 
longer  appears  necessary  to  differentiate 
between  the  parts  of  the  diffusion  flame  in- 
side and  outside  of  an  ammonium  perchlorate 
monopropellant  flame.) 

6)  The  procedure  for  calculation  of  heat  feed- 
back from  the  diffusion  flame  and  the  AP  mono- 
propellant flame  is  modified.  (In  the  BDP 
model,  all  flames  are  considered  to  occur  in 
flame  sheets  at  discrete  distances  from  the 
surface:  in  the  current  model,  the  AP  mono- 
propellant heat  release  is  treated  as  a flame- 
sheet  type  heat  release,  but  the  diffusion 
flame  heat  release  is  considered  to  occur 
over  a finite  range  of  distances  from  the 
propellant  surface.) 

7)  The  distance  (measured  normal  to  the  propel- 
lant surface)  associated  with  oxidizer-binder 
gas  interdiffusion  in  the  presence  of  cross- 
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flow  is  assumed  to  be  reduced  by  a factor, 
sin  8,  where  8 is  the  angle  of  the  resultant 
of  the  crossflow  and  transpiration  velocities 
(at  the  outer  edge  of  the  diffusion  flame 
region)  relative  to  the  surface,  as  in  the 
first  generation  model  described  in  Ref.  1. 

Model  Development 

A major  assumption  made  in  the  BDP  model  (and 
variants  thereof)  is  that  one  may  work  in  terms  of 
an  average  oxidizer-fuel  ratio  for  a given  size  oxi- 
dizer particle.  In  reality,  however,  an  oxidizer 
particle  and  the  fuel  surrounding  it  (and  associated 
with  it)  will  be  receiving  heat  feedback  from  a dif- 
fusion flame  of  strongly  varying  oxidizer/fuel  ratio 
during  its  burning.  As  the  oxidizer  particle  first 
becomes  exposed  to  the  surface,  with  only  its  tip 
showing,  the  local  oxidizer-fuel  ratio  will  be  quite 
low.  As  the  burning  surface  passes  the  equator  of 
the  particle,  however,  the  oxidizer-fuel  ratio  will 
be  comparatively  high,  and  as  the  particle  burns 
out,  the  ratio  will  again  be  low.  Implicit  in  the 
BDP  use  of  an  "average"  oxidizer-surface  planar  in- 
tersectional area  is  the  assumption  that  all  of  the 
highly  non-linear  dependencies  of  burning  rate, 
flame  temperature,  and  consequently  heat  feedback 
from  the  diffusion  flame  can  be  linearly  averaged 
over  the  range  of  the  variations  during  regression 
of  the  propellant  surface  through  the  oxidizer. 
Things  may  work  out  this  way,  but  this  appears  to 
this  author  to  be  a somewhat  risky  a priori  assump- 
tion. Accordingly,  in  this  model  (limited  thus  far 
to  unimodal  oxidizer)  an  attempt  is  made  to  allow 
for  the  variation  in  local  oxidlzer/fuel  ratio 
associated  with  the  burning  of  ' n individual  oxi- 
dizer particle  due  to  the  variation  in  relative  oxi- 
dizer-fuel surface  intersectional  areas  as  the  sur- 
face regresses  through  the  particle. 

In  deciding  how  to  treat  this  variation  (or, 
indeed,  whether  to  treat  it)  one  must  first  address 
the  question  of  propellant  surface  and  subsurface 
response  to  variation  in  heat  feedback  flux  from  the 
varying  oxidlzer/fuel  gas-phase  diffusion  flame.  If 
the  burning  rate  response  is  very  slow,  such  varia- 
tions in  feedback  flux  are  damped  out  and  the  aver- 
aging procedure  of  BDP  is  probably  adequate.  If, 
on  the  other  hand,  response  of  burning  rate  varia- 
tions to  heat  feedback  flux  variations  is  suffic- 
iently fast,  one  may  use  quasi-steady  state  calcu- 
lations of  the  burning  rate  at  each  fuel/oxldlzer 
area  ratio  during  the  regression  of  the  burning  sur- 
face through  the  particle  and  then  properly  average 
these  to  arrive  at  an  average  burning  rate.  In  be- 
tween these  extremes  lies  great  difficulty.  A 
transient  heat  conduction  program  allowing  for  sur- 
face ablation  was  employed  to  examine  the  response 
of  ablation  rate  to  variation  in  heat  flux  to  the 
surface.  Variations  in  heat  flux  up  to  106  cal/cm2 
sec2  (corresponding  to  approximate  doubling  of  heat 
feedback  flux  from  a typical  steady-state  value  in 
0.50  msec,  the  time  required  for  a propellant  burn- 
ing at  2 cm/sec  to  regress  10  microns)  were  examined. 
In  all  cases,  the  burning  rate  response  was  found 
to  track  the  feedback  flux  variation  within  10  per- 
cent. Accordingly,  it  was  concluded  that  use  of  a 
quasl-steady-state  approach  to  calculation  of  pro- 
pellant burning  rate  at  various  oxidlzer/fuel  ratios 
associated  with  different  intersections  of  the  pro- 
pellant burning  surface  with  a given  oxidizer  par- 
ticle would  not  be  seriously  in  error. 


concluded  that  one  can  use  a quasl-steady-state 
approach  to  calculating  burning  rate  as  a function 
of  the  ratio  of  planar  areas  of  oxidizer  and  assoc- 
iated fuel  Intersected  by  the  regressing  surface, 
one  is  next  faced  with  the  question  of  how  to  cal- 
culate the  distribution  of  these  areas.  Since  com- 
posite propellants  are  normally  quite  highly  loaded 
with  solid  oxidizer  in  the  rubber  fuel  binder,  and 
since  with  unimodal  oxidizer  propellants  the  desire 
for  these  high  loadings  tends  to  lead  to  loadings 
approaching  maximum  theoretical  loading,  it  was 
decided  that  as  a reasonable  approximation,  one 
might  assume  a regular  packing  of  oxidizer  crystals 
in  the  binder  corresponding  to  the  arrangement  of  a 
cubic  closest  packing  array,  though  with  the  spac- 
ing larger  than  that  for  a true  cubic  closest  pack- 
ing, corresponding  to  less  than  100  percent  of  theo- 
retical maximum  loading.  Simple  geometrical  con- 
siderations then  permit  one  to  calculate  the  char- 
acteristic lattice  dimension  D^  (where  lattice  8 pac- 
ings in  three  mutually  orthogonal  planes  are  given 
by  Di,  0.866  D^,  and  0.82  Dj)  as: 


B /Ot237\1/3n  /0.737[HBO/pox-V(l-«FO)/pp]pox 
Dl"\  VL0  / o\  WPO 


(1) 


It  is  arbitrarily  assumed  that  the  propellant  burns 
in  the  direction  in  which  the  planes  of  oxidizer  are 
separated  by  0.82  Dj.  This  distance  is  broken  up 
into  equally  spaced  Increments  and  .straightforward 
geometrical  relations  are  then  used  to  calculate  the 
planar  intersection  area  of  the  burning  surface  with 
the  oxidizer  (AF0X)  and  its  associated  fuel  planar 
area  (AFU)  at  each  of  the  intersection  planes,  with 
the  assumption  that  wherever  two  layers  of  oxidizer 
overlap  the  fuel  is  apportioned  between  them  in  the 
ratio  of  their  planar  surface  Intersection  areas. 

The  result  of  these  calculations  is  a table  of  planar 
oxidizer-surface  intersectional  area  (APOX)  and 
associated  fuel  surface  area  (AFU)  versus  distance 
of  the  Intersection  plane  from  the  initial  top  of 
the  particle  (XDTDP) . Results  of  a typical  calcu- 
lation are  presented  in  Table  1.  Burning  rates  for 
each  of  these  conditions  (starting  at  the  top  of 
the  particle  since  one  must  allow  for  different  re- 
gression rates  of  fuel  and  oxidizer)  are  then  cal- 
culated as  described  below  and  an  averaging  pro- 
cedure, also  described  below,  is  then  used  to  cal- 
culate an  average  propellant  burning  rate. 


Next  let  us  address  the  question  of  calculation 
of  propellant  burning  rate  at  each  of  the  conditions 
defined  by  the  various  distances  of  the  burning  sur- 
face Intersection  plane  from  the  top  of  the  oxidizer 
particle,  as  listed  for  the  example  in  Table  1. 
First,  since  as  mentioned  above,  different  oxidizer 
and  fuel  regression  rates  are  to  be  allowed,  one 
must  address  rather  carefully  the  questions  of  sur- 
face geometry  and  mass  conservation  at  the  surface. 
In  this  model,  as  in  the  BDF  model,  the  fuel  is 
assumed  to  regress  in  a planar  manner,  and  the  oxi- 
dlzer-fuel  surface  is  forced  to  be  continuous  at 
their  Intersection.  These  restrictions,  coupled 
with  the  fact  that  the  linear  regression  rates  of 
fuel  and  oxidizer  perpendicular  to  their  directions 
of  regression  are  allowed  to  differ,  force  the  oxi- 
dizer surface  to  assume  a curved  shape  as  it  re- 
gresses. Oxidizer  mass  fluxes  may  be  expressed 
relative  to  either  the  actual  curved  surface  area 
or  the  planar  projection  of  this  area,  the  two 
values  being  related  by: 
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A (APOX)  =A  (AS0X)  = -2*  (AS0X) 

ox,p  ox,s  Pox 
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As  mentioned  earlier,  this  model  is  presently 
limited  to  unimodal  oxidizer  particle  size.  Having 


(2,2a) 


The  area-average  mass  flux  of  fuel  and  oxidizer 
normal  to  the  mean  regression  plane  is  given  by: 

. Afuel<AFU>  + Aox.n<ATOX> 

” AFU  + A POX 

Afuel(AFU)  + Aox,.(ASOX) 

AFU  + A POX 


(3) 


It  is  important  to  know  the  value  of  ASOX  at  each 
plane  since  Che  Arrenhius  expression  relating  oxi- 
dizer mass  flux  to  surface  temperature  must  be  writ- 
ten in  terms  of  AoX,s  t0  be  meaningful: 


A 


*ox,s  = Box  exP  (-Eox/RV  <4) 

similar  expression  for  the  fuel  pyrolysis  rate: 

“fuel  = Bfuel  exp(-Efuel/RV  (5> 


enables  one  to  calculate  the  ratio  of  oxidizer  and 
fuel  regression  rates  as  a function  of  surface  tem- 
perature : 

rox  _ mox.s  Pfuel 

rfuel  *fuel  Box  ... 

lb; 

p,  , B exp(-E  /RT  ) 

^fuel  ox  r ox  s 

= p B.  ,exp(-E,  , /RT  ) 

'ox  fuel  r fuel  s 

There  is  considerable  uncertainty  as  to  best  values 
to  be  used  for  Box,  Bfuei,  Eox,  and  Efuel*  thus, 
parametric  study  of  the  effects  of  these  values  is 
required.  Note  that  it  has  been  assumed  here  that 
the  oxidizer  and  fuel  surface  temperatures  are  equal. 
This  is  probably  not  a particularly  good  assumption, 
but  relaxing  it  requires  a rather  complex  three- 
dimensional  heat  transfer  analysis. 


Now,  how  does  one  go  about  calculating  ASOX  for 
succeeding  regression  intervals  through  the  oxidizer 
particle?  First,  it  is  assumed  (approximated)  that 
at  the  first  increment  after  the  tip  of  the  particle 
becomes  exposed  (in  Table  1,  when  the  distance  from 
the  top  of  the  particle  is  0.119  microns)  the  oxi- 
dizer surface  is  planar.  The  procedure  outlined 
below  for  calculation  of  burning  rate,  given  the 
local  oxidlzer/fuel  area  ratio,  is  then  used  to  cal- 
culate the  oxidizer  and  fuel  linear  regression  rates 
under  the  conditions  given  for  this  first  Increment. 
The  fuel  regression  rate  is  then  used  to  calculate 
the  time  for  the  regressing  fuel  to  reach  the  second 
increment  (distance  from  the  initial  particle  top 
of  1.017  microns  in  Table  1)  as: 

TAU2=  (XDT0P2  - XDTOPj)  /r£  t (7) 

The  distance  which  the  center  of  the  AP  particle 
peak  regresses  in  that  time  is  then  calculated  as: 

A(OEIiOX)  = rQx(TAU2)  (8) 

Similar  procedures  are  followed  for  each  succeeding 
Increment,  yielding  for  each  XDTOP  (distance  of  fuel 
surface  from  the  initial  top  of  the  oxidizer  par- 
ticle) a value  of  DEIOX  (distance  of  the  center  of 
the  oxidizer  surface  from  the  initial  top).  The 
geometrical  method  depicted  in  Figure  1 (and  dis- 
cussed in  detail  in  References  3 and  6)  then  permits 
calculation  of  the  actual  curved  oxidizer  surface 
area,  ASOX,  as: 

ASOX  = TT[ (XDTOP) Do  - 2 (XDTOP)  (DEIOX)  + DEIOX2]  (9) 

For  the  calculation  of  the  length  of  the  col- 
umnar diffusion  flame  between  fuel  gases  from  the 


binder  pyrolysis  and  oxidizer  gases  from  the  oxi- 
dizer decomposition  (discussed  later)  several  other 
parameters  associated  with  the  surface  configuration 
of  the  oxidizer  particle-associated  fuel  combination 
at  each  Increment  must  be  calculated.  First,  the 
combined  radius  of  the  oxidizer  and  binder  gas 
streams  (in  this  model,  a modified  Burke-Schumann 
analysis  with  a fuel  annulus  surrounding  an  oxidizer 
gas  core  is  employed  for  the  columnar  diffusion 
flame  calculation)  is  calculated  as: 

Rjg  = ./(AFU  + APOX)  /TT  (10) 

In  line  with  the  requirement  in  the  Burke-Schumann 
analysis  that  the  linear  velocity  of  the  fuel  and 
oxidizer  streams  have  the  same  Initial  value,  the 
oxidizer  and  fuel  gas  streams  leaving  the  surface 
are  assused  to  adjust  their  areas  quickly  from  the 
planar  area:,  of  the  solids  to  meet  this  requirement. 
Since  the  temperature  at  the  surface  is  assumed  to 
be  the  same  for  the  fuel  and  oxidizer,  and  pressures 
are  equal,  this  leads  to  an  expression  for  the 
radius  of  the  inner  oxidizer  gas  jet  of: 


/( 


1 + 


d,tuei(AFU)(HW)ox  \ 
Aox,s(AS0X)(MW)fuel/ 


(U) 


The  linear  gas  velocity  away  from  the  surface, 
also  required  in  the  modified  Burke-Schumann  analy- 
sis, as  well  as  for  calculation  of  characteristic 
reaction  distances  (products  of  reaction  times  and 
this  velocity)  is  calculated  as: 


*nx.s(AS°X)RTsurf 

<™w2p 


(12) 


Finally,  the  ratio  of  the  molar  fuel/oxidlzer 
ratio  to  stoichiometric  molar  fuel/oxidizer  ratio 
for  the  combined  fuel  and  oxidizer  streams  ($)»  also 
required  in  the  modified  Burke-Schumann  analysis, 
is  calculated  as: 


t = 


(13) 


Since  the  initial  pressure  and  temperature  of  the 
fuel  and  gas  streams  are  the  same,  the  concentration 
ratio  C2/C1  may  be  replaced  by  a ratio  of  mole  frac- 
tions Yj>/Yl.  With  subsurface  reactions,  the  fuel 
mole  fraction  is  reduced  by  a factor  (1-R)  from  its 
value  in  the  absence  of  any  subsurface  reactions, 
while  the  oxidizer  mole  fraction  is  reduced  from  its 
no-subsurface-reaction  value  by  a factor  (l-o) . With 
these  substitutions,  Equation  13  becomes: 


i = 


Y2.o(1  - B)(RBS2  ~ LBS2> 

1 Yi,o(1  - •>  h»s2 


(14) 


The  burning  rate  of  the  propellant  at  any  given 
set  of  oxidlzer/fuel  conditions  (any  regression  in- 
crement) is  controlled  by  heat  released  (exothermic 
reactions)  at  various  locations.  In  this  model,  we 
consider  three  principal  heat  release  zones:  (1) 
heat  release  In  a thin  subsurface  zone  quite  near 
(and  including)  the  propellant  surface;  (2)  heat 
release  In  the  gas-phase  above  the  propellant  from 
ammonium  perchlorate  decomposition  products  burning 
as  a monopropellant;  and  (3)  heat  release  from  a 
diffusion  flame  between  AP  decomposition  (and  mono- 
propellant flame)  products  and  fuel  vapor  released 


by  binder  pyrolysis 


The  subsurface /surface  heat  release  is  calcu- 
lated by  an  iterative  process,  coupled  with  the  re- 
mainder of  the  model,  in  which  an  estimate  of  the 
subsurface  temperature  profile  Is  made  and  substi- 
tuted into  an  Arrenhius  rate  expression  represent- 
ing subsurface  heat  release  rate  data  measured  by 
Waesche  and  Wenograd^,  which  is  then  integrated 
from  the  surface  deep  into  the  unburned  propellant 
to  obtain  the  total  subsurface  heat  release  per 
unit  mass  of  propellant.  This  procedure  differs 
markedly  from  that  of  the  BDP  model,  in  which  the 
amount  of  subsurface  heat  release  per  unit  mass  of 
propellant  is  assumed  to  be  a constant.  Independent 
of  such  parameters  as  burning  rate,  and  is  included 
with  the  binder  heat  of  vaporization.  Since  the 
subsurface  temperature  profile  steepens  rapidly  with 
increasing  burning  rate,  while  surface  temperature 
increases  with  burning  rate,  our  procedure  results 
in  the  subsurface  heat  release  per  unit  mass  of  pro- 
pellant varying  with  burning  rate.  As  will  be  dis- 
cussed later,  the  surface  energy  balance  in  this 
model  is  written  with  the  surface  area  of  the  oxi- 
dizer and  associated  fuel  as  the  basis:  thus  all 
terms  appear  in  the  units  of  energy/time.  For  book- 
keeping convenience,  the  surface/subsurface  heat 
release  term  is  written  as: 


4sub 
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It  is  assumed  that  a stoichiometric  amount  of  fuel 
is  reacted  with  the  oxidizer  in  these  surface/sub- 
surface reactions.  Thus  the  fraction  of  fuel  re- 
acted in  these  reactions,  8,  is  given  by: 

m (ASOX) 

e = SMRBO  or  (16) 
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Based  upon  Waesche  and  Wenograd  subsurface  reaction 
rate  data,  the  fraction  of  oxidizer  reacted  per  unit 
time  is  given  for  AP/CTPB  systems  as  a function  of 
temperature  by: 


R = B.  . exp(-E  . /RT) 
a Sub  r sub 


(17) 


with  preliminary  examination  of  their  data  indicat- 
ing that  the  activation  energy  is  approximately 
40,000  calories /mole,  while  the  pre-exponential  is 
approximately  2.5  (10") . (As  will  be  discussed 
later,  however,  use  of  these  values  leads  to  very 
little  predicted  subsurface  heat  release  which  in 
turn  leads  to  underprediction  of  burning  rate  for 
formulations  containing  large  oxidizer  particles: 
accordingly,  somewhat  higher  values  of  the  pre-ex- 
ponential factor  were  also  examined.)  The  unper- 
turbed (uncoupled)  subsurface  temperature  profile 
is  given  by: 


T = (T 
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Substitution  of  Equation  18  into  Equation  17  and 
Integration  of: 
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(19) 


where  x'  represents  a distance  below  the  surface  at 
which  the  reaction  rate  has  dropped  to  1 percent  of 
its  surface  value,  yields,  with  use  of  several 
approximations : 6 
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Thus  Equation  20  relates  the  mass  fraction  of  oxidi- 
zer reacting  exothermically  at  or  below  the  surface 
to  the  surface  temperature  and  linear  regression 
rate  of  the  oxidizer. 


As  regards  gas-phase  heat  release  zones,  a two- 
flame  approach  was  chosen  for  this  model,  the  two 
flames  being  an  AP  monopropellant  flame  and  a col- 
umnar diffusion  (Burke-Schumann)  flame.  The  reasons 
that  a two-flame  rather  than  a three-flame  model 
(as  the  BDP)  was  chosen  were: 


1)  Mathematical  simplication. 


2) 


Lack  of  apparent  difference  in  a diffusion 
flame  between  AP  decomposition  products  and 
fuel  and  a flame  between  AP  monopropellant 
flame  products  and  fuel.  In  both  cases,  the 
overall  stoichiometry  Is  the  same  Iftnce, 
while  AP  decomposition  products  bring  more 
oxidizer  into  a binder  fuel  stream  than  do 
AP  monopropellant  flame  products , they  also 
bring  more  fuel,  with  the  result  that  the 
overall  mixture  ratio  at  a given  point  is 
nearly  the  same. 


3) 


Provisions  were  made  in  calculation  of  the 
AP  monopropellant  heat  release  to  allow  for 
consumption  of  reactants  for  that  flame  in 
part  of  the  columnar  diffusion  flame  which 
occurred  inside  the  AP  flame. 


Three  distance  parameters  are  important  in  cal- 
culating heat  feedback  from  these  gas  flames  to  the 
propellant  surface.  These  are  pictured  in  Figure  2. 
These  distances  are  designated  as  FH90  sin  8,  L^p 
and  Lrx-  FH90  refers  to  the  distance  associated 
with  completion  of  90  percent  of  the  mixing  of  fuel 
and  oxidizer  gas  products.  (If  there  were  no  re- 
action delay,  this  would  be  equivalent  to  the  90 
percent  heat  release  point.)  L]{x  refers  to  the  re- 
action distance  associated  with  the  binder  gas-oxi- 
dizer gas  flame,  and  Lap  refers  to  the  reaction  dis- 
tance associated  with  the  monopropellant  AP  product 
flame  (both  being  characteristic  reaction  times 
multiplied  by  the  gas  velocity  away  from  the  sur- 
face) . As  discussed  in  Reference  1,  flame  bending 
associated  with  crossflow  is  assumed  to  reduce  the 
distance  from  the  surface  to  the  end  of  the  columnar 
diffusion  heat  release  (90  percent  point)  by  reduc- 
ing FH90  to  FH90  sin  9,  measured  perpendicular  to 
the  surface,  where  0 is  the  angle  between  the  sur- 
face and  the  resultant  vector  of  the  transpiration 
and  crossflow  velocities  at  the  outer  edge  of  dif- 
fusion flame  zone.  FH90  is  calculated  as  a function 
of  various  parameters  using  a modified  Burke-Schumann 
analysis  as  described  below.  A series  of  calcula- 
tions of  FH90  as  a function  of  these  parameters  was 
performed  externally  and  correlations  of  the  results 
were  used  in  the  final  program.  In  this  model,  it 
is  assumed  that  the  fraction  of  planar  projection 
of  surface,  AP0X/ (AP0X  + AFU)  receives  flux  from 
both  AP  monopropellant  and  columnar  diffusion  flames 
(the  latter  at  the  adiabatic  flame  temperature,  Tf, 
which  is  a function  of  the  oxidlzer/fuel  ratio, 
Aox,s[AS0X]/i5ifUel [AFU ])  while  the  remaining  fraction 
of  the  surface  receives  flux  only  from  this  diffus- 
ion flame;  however,  these  fluxes  are  assumed  to 
smear  out  uniformly  In  the  propellant.  Thus,  the 
total  heat  flow  from  the  gas -phase  heat  release 
zones  is  given  as: 
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Heat  release  from  the  AP  monopropellant  flame  la 
assumed  to  be  a planar  one,  resulting  In  a discon- 
tinuity in  the  temperature  derivative  at  its  point 
of  release,  while  the  columnar  diffusion  flame  is 
assumed  to  release  its  heat  uniformly  between  x - 
Lrx  and  x - Lpx  + FH90  sin  9.  (Actually,  detailed 
examination  of  the  Burke -Schumann  flame  structure 
indicates  that  this  heat  release  should  be  more 
heavily  weighted  toward  the  inner  side  of  this  zone 
- this  will  be  discussed  further  later  in  the  paper.) 
With  use  of  multitudinous  algebraic  manipulations, 
we  arrive  at  the  following  expressions  for  the  heat 
fluxes  at  the  surface  (allowing  for  reduction  of 
reactants  available  for  the  AP  monopropellant  re- 
action by  occurrence  of  some  diffusion  flame  react- 
ion closer  to  the  surface) : 
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(z2e  L-  e + 1) 

Z1  = “V^RX  + 0,90  Sln  ®)/X 

z,  = AC  FH90  sin  0/A 
2 P 

z3  = ACp(FH90  sin  0 + L^-  LAp)/A 

Next,  let  us  consider  the  calculation  of  the 
distances  l^x,  lAP,  and  FH90  (and  FH90  sin  9).  The 
distance  FH90,  which  is  calculated  from  a modified 
Burke-Schumann  columnar  diffusion  flame  analysis 
(modified  to  allow  for  axial  diffusion)  is  defined 
as  the  distance  from  the  starting  plane  at  which  90 
percent  of  the  fuel  (for  oxidizer-rich  cases)  or  90 
percent  of  the  oxidizer  (for  fuel-rich  cases)  will 
be  consumed,  assuming  infinite  reaction  kinetics. 
This  definition  of  the  characteristic  diffusion 
distance  differs  from  that  of  the  BDP  model  where 
the  characteristic  distance  is  defined  as  the  dis- 
tance from  the  starting  plane  to  the  point  of  clo- 
sure of  the  flame  over  the  oxidizer  (fuel-rich 
cases)  or  the  fuel  (oxidizer-rich  cases) , One  ser- 
ious problem  with  use  of  the  flame  closure  point 
to  define  the  characteristic  distance  is  that  it 
has  a singularity  for  stoichiometric  situations: 
that  is,  for  stoichiometric  oxidizer-fuel  ratio, 
the  flame  does  not  close  and  this  characteristic 
distance  goes  to  infinity.  Since  most  of  the  heat 
is  still  released  fairly  close  to  the  surface,  this 
latter  definition  of  a characteristic  diffusion  dis- 
tance leads  to  seriously  misleading  results  as  re- 
gards heat  feedback  to  the  propellant  surface  at 
near  stoichiometric  oxidizer-fuel  ratios:  at  oxi- 
dizer-fuel ratios  far  from  stoichiometric,  FH90 
and  the  distances  associated  with  flame  closure 


differ  only  slightly.  Details  of  the  calculation 
of  FH90  are  discussed  in  Reference  6.  This  para- 
meter is  a function  of  four  independent  parameters: 
the  oxidizer  jet  radius,  the  equivalence  ratio,  the 
ratio  of  the  outer  radius  of  the  fuel  annulus  to 
the  oxidizer  jet  radius,  and  the  ratio  of  dlffusi- 
vlty  to  transpiration  velocity  (evaluated  at  T 
(Ts+Tf)/2). 

FH90  - fjd*.  D/Vgas,  4,  Rbs/Lbs)  (24) 

An  extensive  set  of  calculations  covering  a wide 
range  of  each  of  these  variables  was  carried  out 
externally  to  the  final  computer  code  and  tabula- 
tions and  correlations  of  the  results  were  built 
into  the  final  combustion  model. 

The  reaction  distances,  Lgx  and  LAp>  are  cal- 
culated as  the  products  of  reaction  times  and  gas 
velocity  away  from  the  surface.  Using  the  Zeldovlch 
approach  for  premixed  flame  analysis  along  with  sev- 
eral minor  approximations  which  will  not  be  detailed 
here,  we  arrive  at  for  the  oxidizer-fuel  reaction 
distance,  Lgy: 

, K0FV»s.,urf<1+»>2Tf2<eZ^VT.QF^F^ 
hix-  p(n-l)t  (25) 

For  the  anmonium  perchlorate  reaction  distance, 
we  find,  neglecting  variation  of  temperature  at  the 
AP  heat  release  site  (probably  not  a very  bad  approx- 
imation due  to  the  low  activation  energy  associated 
with  the  airnnonia-perchloric  acid  reaction): 
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The  same  approach  to  calculation  of  sin  9 as 
used  in  the  previously  mentioned  first  generation 
erosive  burning  model*  was  used  in  this  model,  the 
resulting  equations  being: 

Tf_ 

Ugas,x=FH90  sin  0 Vgas,surf  Tg  ^273 

Y+  IO-5(FH9°  slnS)  ^crossflow)0’ V,9p0,9 

D0*1  ,(tT  ♦ Tf]/2)1’8  ( > 

channel  s f 

(FH90  in  microns;  U in  ft/sec;  D in  ft;  T in  °K; 

P in  atm) 

U+  = Y+  for  Y+  <5 

U+  = -3.05  + 5.00  lnY+  for  5 < Y+  < 30  (29) 

U+  = 5.5  + 2.5  lnY+  for  Y+  >30 


crossflow, x = FH90  sin  9 
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gln  e = gas.x  = FH90  sin  9 (31) 

+ u2 

’ gas,x=FH90  sin  0 crossf low,x=FH90  sin  9 

As  mentioned  earlier,  the  final  flame  tempera- 
ture, Tf , depends  on  the  relative  flow  rates  of  fuel 


and  oxidizer  gases  at  each  calculatlonal  increment 
during  progression  of  the  propellant  surface  through 
the  oxidizer.  Accordingly,  a table  of  flame  temper* 
ature  calculated  with  a thermochemical  equilibrium 
program  versus  a parameter  representing  the  relative 
flow  rates  is  generated  for  the  propellant  system  of 
interest  (e.g.,  HTPB-AP)  and  included  as  a tabular 
look-up  in  the  final  program  in  the  form: 
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In  addition,  the  product  gas  heat  capacity  is 
somewhat  dependent  upon  this  parameter  and  an  addi- 
tional tabular  look-up,  based  on  thermochemical  cal- 
culations, is  Included. 


At  this  point,  we  have  26  equations  (2, 2a, 3-6, 
10-12,  14-16,20-33)  in  27  unknowns  (&oX|p,  Aq*  s, 
rox»  mfuel»  Ts,  rfuei,  RbS»  ^gas,surf*  Cpi  n, 
Pt  9sub>  'gas*  ^series  flames*  ^diffusion  flame*  ^BS, 
IAP,  Lrx,  FH90,  »,  Tf,  Ugas,x=FH90  sin  6>  T4" , IT*-, 
ucrossflow.x=FH90  sin  e) • for  closure  of  the  pro- 
blem, we  finally  write  an  energy  balance  at  the 
propellant  surface  as: 
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The  resulting  27  equations  (some  of  which,  as 
mentioned,  are  tabular  look-ups  or  correlations  of 
one  parameter  as  a function  of  others)  are  solved 
simultaneously  in  a computer  program  for  each  given 
set  of  values  of  ASOX  and  AFU  associated  with  each 
surface  regression  Increment.  Among  the  outputs  of 
each  solution  are  values  for  rox,  and  rfuei  which 
are  used  in  calculation  of  ASOX  from  the  known  AFOX 
for  the  succeeding  Increment  via  Equations  7-9. 

In  the  solution  procedure,  the  thermal  conduct- 
ivity of  the  gas  and  the  ratio  of  dlffuslvlty  to  gas 
velocity  (one  of  the  Independent  parameters  in  Equa- 
tion 24) , both  proportional  to  the  square  root  of 
temperature,  are  evaluated  at  the  average  of  the 
flame  and  surface  temperature. 

As  the  program  is  stepped  through  the  succeed- 
ing Increments  of  fuel  plane  distance  from  the  ini- 
tial top  of  the  oxidizer  particle  (see  Table  1) , the 
oxidizer  surface  will  either  assume  a protruding 
bulge  or  a depression  relative  to  the  planar  fuel 
around  it,  depending  upon  whether  the  oxidizer  lin- 
ear regression  rate  is  slower  or  faster  than  the 
binder  regression  rate.  This  raises  interesting 
questions  regarding  the  "end-game"  if  the  particle 
burns  out  before  the  surrounding  fuel.  (Geometrical 
considerations  show  that  the  inverse  problem  cannot 
occur  as  long  as  increment  sizes  are  kept  suffici- 
ently small.)  In  this  case,  there  is  no  oxidizer 
to  bum  with  the  surrounding  fuel  in  succeeding  in- 
crements and  the  burning  rate  is  set  equal  to  zero 
for  these  remaining  Increments.  Three  different 
approaches  have  been  taken  to  calculating  the  aver- 
age propellant  burning  rate  from  the  information 
obtained  during  the  procedure  of  stepping  through 


the  increments  of  regression  of  the  fuel  planar  sur- 
face past  the  oxidizer  particle  (XDTOP  increments). 
In  the  first  of  these,  the  burning  rate  is  calcu- 
lated by  statistically  averaging  all  of  the  oxidizer 
mass  fluxes  and  fuel  mass  fluxes  over  the  Increments 
and  dividing  by  the  propellant  density: 

, ox.p.j  j fuel,)  j 
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while  in  the  second  approach  it  is  calculated  by 
statistically  averaging  all  of  the  oxidizer  mass 
fluxes  and  then  dividing  by  the  overall  oxidizer 
mass  fraction  and  the  propellant  density: 

r“*’  Wu-t™  "V 

The  fact  that  these  two  procedures  do  not  always 
give  the  same  result  (though  the  differences  are 
generally  small)  is  tied  in  with  the  "end-game" 
problem  mentioned  above.  If  the  oxidizer  burns  out 
before  the  fuel  plane  reaches  the  bottom  of  the 
oxidizer,  mass  fluxes  for  succeeding  Increments  are 
set  equal  to  zero  in  the  procedure  currently  used. 

Not  only  does  this  result  in  different  answers  by 
the  two  above  procedures,  but  it  also  pulls  the 
average  burn  rates  down.  One's  first  temptation  is 
to  simply  perform  the  sunmlng  procedure  over  just 
these  increments  for  which  a burning  rate  is  calcu- 
lated, both  in  the  numerator  and  denominator  of 
Equations  35  and  36,  but  it  is  not  clear  whether 
or  not  this  is  more  physically  realistic  than  the 
procedure  of  summing  over  all  increments,  with  burn- 
ing rate  set  equal  to  zero  for  increments  in  which 
there  is  no  oxidizer  for  the  fuel.  A third  pro- 
cedure of  calculating  average  burning  rate  was  dev- 
eloped which  basically  does  take  this  alternate 
approach,  however,  though  in  a slightly  different 
manner.  In  this  procedure,  the  burning  rate  is  cal- 
culated by  dividing  the  oxidizer  particle  diameter 
by  the  sum  of  the  times  required  for  each  increment 
until  the  bottom  of  the  oxidizer  particle  is  reached: 
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This  procedure  begs  the  question  of  what  happens  to 
the  fuel  "left  over"  when  the  oxidizer  particle 
bums  out  before  the  fuel.  Physically  we  can  per- 
haps just  assume  that  it  somehow  flakes  off.  This 
question  needs  to  be  addressed  further.  This  third 
approach  does  allow  for  the  fact  that  the  particle 
will  spend  more  time  at  regression  Increments  where 
the  burning  rate  is  lower,  while  the  first  two  pro- 
cedures involve  an  implicit  assumption  that  each 
of  the  regression  Increments  is  equally  likely. 


During  the  past  year,  Atlantic  Research  has 
collected  extensive  experimental  burning  rate-pres- 
sure data  in  the  presence  of  and  absence  of  cross- 
flow  for  three  non-metalllzed  composite  formulat- 
ions containing  unimodal  acirvonium  perchlorate  (AP) 
oxidizer.  All  three  of  these  formulations  consist 
of  73  velght  percent  AP  and  27  weight  percent  hydro- 
xytermlnated-polybutadiene  (HTPB)  binder  (with  a 
trace  of  carbon  black  to  opacify  them) . One  of  the 
formulations,  designated  as  Formulation  4685,  con- 
tains 5 micron  diameter  AP;  the  second,  Formula- 
tion 4525,  contains  20  micron  diameter  AP;  and  the 
third.  Formulation  5051,  contains  200  micron  diameter 


AP.  Preliminary  testing  of  the  model  described 
above  lias  been  carried  out  against  burning  rate  data 
obtained  for  these  three  formulations.  As  may  be 
seen  from  the  above  equations,  numerous  physical 
constants  (c.g.  heat  capacities,  densities,  thermal 
conductivities,  stoichiometry  parameters,  reaction 
heats,  and  rate  constants)  must  be  estimated  for  use 
of  the  model.  Values  for  all  but  four  of  these  par- 
ameters, estimated  for  the  case  of  HTPB  binder  and 
AP  oxidizer,  appear  in  Table  2.  Discussion  of  the 
estimation  of  these  parameters  appears  in  Reference 
6.  Of  these  parameters,  the  coefficient  for  the 
diffusivity  expression  is  the  least  certain,  since 
this  value  strongly  depends  on  the  nature  of  the 
reacting  gas  species  and  their  products.  It  is  esti- 
mated that  this  coefficient  should  lie  somewhere 
between  1.0(10)"®  and  2.5(10)"®:  the  former  value 
was  chosen  for  the  calculations  presented  in  this 
paper.  Four  parameter  values  are  not  listed  in 
Table  2:  these  are  the  gas  reaction  distance  con- 
stants K^p  and  Kqjt,  the  pre-exponential  term  for  the 
subsurface  reaction  Bsub,  and  the  gas-phase  reaction 
order  n.  During  the  course  of  this  study  gas -phase 
reaction  orders  of  1.8  and  2.0  were  used:  at  this 
stage  the  former  value  appears  to  give  better  re- 
sults and  was  used  in  all  the  predictions  presented 
in  this  paper.  Based  on  crude  estimation  procedures 
discussed  in  Reference  6,  it  was  felt  that  K^P  should 
probably  lie  between  approximately  .00003  and  .00015 
atm0-®  sec  for  a reaction  order  of  1.8,  or  .00005 
and  .0003  atm  sec  for  a reaction  order  of  2.0.  A 
very  crude  estimate  of  Kgp  indicates  that  it  should 
lie  between  0.3(10)"!®  and  1.5(10)-13  atm0-®  sec/°K2 
for  a reaction  order  of  1.8.  Various  values  of 
these  two  parameters  within  these  ranges  were  tried 
in  attempts  to  fit  the  experimental  burning  rate 
data  with  one  set  of  constants  over  the  wide  ranges 
of  pressure,  crossflow  velocity,  and  particle  size 
experimentally  studied. 

As  mentioned  earlier,  our  preliminary  examina- 
tion of  Wacsche  and  Wenograd4  subsurface  reaction 
rate  data  indicated  that  the  fraction  of  oxidizer 
reacted  per  unit  time  in  AP/CTPB  (and  presumably 
AP/HTPB)  systems  could  be  expressed  as  Ra  = 2.5(10^® 
g-AOOOO/R'r  Accordingly,  in  our  first  attempts  at 
using  the  model  described  above,  we  employed  a value 
for  Bauj,  of  2,5(10)!!.  With  this  value,  we  were 
able  to  find  a set  of  values  of  Kap  and  K^p  which 
gave  good  agreement  between  theory  and  experiment 
for  zero  crossflow  for  the  5 micron  and  20  micron 
AP  formulations,  as  shown  in  Figures  3 and  4.  How- 
ever, the  model  did  tend  to  underpredict  the  effect 
of  crossflow  and,  more  important,  resulted  in  gross 
underprediction  of  the  burning  rate  and  the  burning 
rate-pressure  slope  for  the  no-crossflow  case  for 
the  200  micron  AP  formulation.  (Figure  5)  In  all 
of  these  cases,  the  fraction  of  oxidizer  and  fuel 
consumed  in  subsurface  reactions  was  predicted  to 
be  quite  small  (on  the  order  of  1 to  2 percent) . 

Study  of  this  problem  of  underprediction  of 
burning  rate  of  the  200  micron  AP  formulation  indi- 
cated that  it  might  be  corrected  by  one  or  both  of 
two  possible  modifications  to  the  model,  only  one 
of  which  has  been  examined  thus  far.  As  mentioned 
earlier,  the  fuel  gas  - oxidizer  gas  diffusion 
flame  is  assumed  to  release  its  heat  uniformly  be- 
tween X = Lgx  and  X - Lgx  + FH90  sin  0.  (Figure  2). 
Detailed  examination  of  reactant  concentration  pro- 
files in  the  diffusion  flame  as  predicted  by  the 
Burke -Schianann  analysis  indicate  that  in  reality 
the  heat  release  should  be  weighted  fairly  heavily 
toward  the  X = Lgx  side  of  this  zone,  tapering  off 


toward  X Lrx+  FH90  sin  0.  To  allow  for  this  in 
the  model,  alternate  expressions  to  Equations  22 
and  23  relating  gas-phase  heat  feedback  fluxes  to 
FH90  sin  0,  Lrx,  and  Lap  would  have  to  be  derived 
and  substituted:  this  has  not  yet  been  done.  It 
can  be  logically  deduced  that  such  a modification 
would  reduce  the  predicted  dependency  of  burning 
rate  on  oxidizer  particle  size,  thus  perhaps  per- 
mitting good  fitting  of  the  5,  20,  and  200  micron 
AP  formulation  data  with  one  set  of  constants. 

The  second  approach  which  should  also  lead  to 
decreased  dependency  of  predicted  burning  rate  on 
oxidizer  particle  size  is  to  increase  the  predicted 
subsurface  heat  release  by  altering  one  or  both  of 
the  constants  in  Equation  17.  Since  the  Waesche 
and  Wenograd^  data  were  obtained  at  relatively  low 
temperatures  (on  the  order  of  600-700°K  compared 
to  typical  predicted  burning  propellant  surface 
temperatures  of  900-1000<>K)  use  of  the  expression 
fitted  to  their  data  does  involve  a dangerous  degree 
of  extrapolation.  Moreover,  while  an  activation 
energy  of  40,000  calories/mole  was  employed  in  mak- 
ing this  extrapolation,  values  ranging  from  40,000 
to  60,000  calories/mole  appear  in  various  reports 
by  these  authors.  If  one  fixes  the  rate  of  reaction 
measured  at  600°K  and  uses  activation  energies  of 
40,000.  50,000,  and  60,000  calories/mole  to  extrap- 
olate to  950°K,  one  predicts  rates  at  that  temper- 
ature in  the  ratio  of  1 to  23  to  480.  Accordingly, 
the  value  of  Bsuj,  was  allowed  to  be  a free  parameter 
in  further  optimization  of  constants  for  the  model. 
Though  the  best  procedure  would  have  been  to  simul- 
taneously alter  Esuj,  so  as  to  give  the  same  predicted 
rate  at  600-6  50°K,  in  this  preliminary  examination 
Esub  was  held  constant  at  40,000  calories/mole. 

Kith  a value  for  Bgub  of  5.2(10)12  (approximately 
20  times  the  value  initially  used)  the  results  shown 
in  Figures  6-8  were  obtained.  As  may  be  seen, 
the  agreement  between  experiment  and  prediction  is 
excellent  for  the  no-crossflow  cases  across  the 
entire  range  of  pressure  and  ptrtlcle  size  studied, 
though  the  effect  of  crossflow  is  still  underpre- 
dicted. With  this  value  of  Baub,  on  the  order  of 
25  to  50  percent  of  the  oxidizer  is  predicted  to  be 
consumed  in  subsurface/surface  reactions. 

It  should  be  emphasized  that  these  results  are 
quite  preliminary.  Further  examination  of  the  sub- 
surface reaction  rate  term  and  of  the  distribution 
of  heat  release  in  the  fuel-oxidizer  diffusion  flame 
will  be  carried  out.  In  addition,  the  question  of 
choice  of  a proper  averaging  procedure  for  bum  rate 
versus  regression  increment  will  be  further  examined. 


A model  for  prediction  of  burning  rate-pressure- 
crossflow  velocity  relationships  for  non-metallized 
composite  propellants  containing  unimodal  oxidizer, 
given  only  composition  and  particle  size,  has  been 
developed.  Data  for  three  such  propellants  of 
identical  composition  (73/27  AP/HTPB)  containing 
different  sizes  of  unimodal  oxidizer  (5,  20,  and 
200  micron  diameter)  have  been  used  for  optimization 
of  three  free  constants  in  the  model.  In  a first 
effort,  only  two  constants,  those  relating  gas  flame 
reaction  distances  to  gas  velocity  away  from  the 
surface  and  pressure,  were  varied  in  a search  for 
optimum  values,  the  reaction  constant  for  exothermic 
subsurface  reactions  being  held  at  a value  based  on 
extrapolation  of  DSC  data  obtained  by  Waesche  and 
Wenograd1^.  In  this  case  a set  of  constants  was 
found  which  yielded  good  agreement  between  experiment 
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and  theory  for  the  no-crossflow  case  for  the  5 and 
20  micron  diameter  AP  formulations , but  seriously 
underpredicted  burning  rates  for  the  200  micron  AP 
formulation.  Examination  of  the  problem  indicated 
that  either  use  of  a higher  rate  constant  for  the 
subsurface/surface  exothermic  reactions  or  a more 
accurate  description  of  heat  release  distribution 
in  the  fuel/oxidizer  diffusion  flame  above  the  pro- 
pellant surface  should  tend  to  reduce  the  sensiti- 
vity of  predicted  burning  rate  to  oxidizer  particle 
size,  as  desired.  Thus  far,  only  the  former  ap- 
proach has  been  examined  and  has  been  found  to  per- 
mit selection  of  a set  of  constants  which  gives  good 
agreement  between  predicted  and  measured  burning 
rate  versus  pressure  curves  for  all  three  particle 
sizes  for  the  case  of  no  crossflow.  In  all  cases, 
however,  the  predicted  sensitivity  of  burning  rate 
to  corss>flow  velocity  is  somewhat  less  than  observed. 
It  must  be  stressed  that  these  results  are  very  pre- 
liminary and  that  further  examination  and  optimiza- 
tion of  the  model  is  in  progress. 

Nomenclature 


AFU 

APOX 

ASOX 

Bfuel 

Box 

Bsub 

CP 

Cp.fuel 

cp,ox 


D 

Dl 

Do 

^channel 

DELOX 

EACT.OF 

Efuel 

Box 

Esub 

FH90 

1 

kAP 

Kof 

lap 

EbS 

LrX 

A 


Fuel  surface  area 

planar  projection  of  exposed  oxidizer 

particle  surface  area 

Total  curved  oxidizer  exposed  surface 

area 

Pre-exponential  in  fuel  pyrolysis  mass 
flux  expression 

Pre-exponential  in  oxidizer  sublimation 
mass  flux  expression 

Pre -exponential  for  subsurface  reaction 
rate  equation 

Gas  heat  capacity  (function  of  oxidizer/ 
fuel  ratio) 

Solid  fuel  heat  capacity 
Solid  oxidizer  heat  capacity 
initial  net  (excess  of  oxidizer  over 
fuel  molecules)  molar  concentration  of 
oxidizer  in  the  oxidizer  decomposition 
product  gases 

initial  molar  concentration  of  fuel  in 
the  binder  pyrolysis  product  gases. 

Gas  diffusivity 

Characteristic  lattice  dimension  in  reg- 
ular oxidizer-particle  array 
Oxidizer  particle  diameter 
Flow  port  hydraulic  diameter 
Distance  of  center  of  oxidizer  surface 
from  initial  oxidizer  peak  (See  Fig.  1) 
Activation  Energy  for  the  fuel-oxidizer 
gas  reaction 

Activation  Energy  for  fuel  pyrolysis 
Activation  Energy  for  oxidizer  sublima- 
tion 

Activation  Energy  for  subsurface  reaction 
rate  equation 

Distance  required  for  mixing  of  fuel  and 
oxidizer  gas  streams 
Stoichiometric  moles  of  fuel  per  mole 
of  oxidizer 

Constant  in  expression  for  oxidizer  mono- 
propellant reaction  distance 
Constant  in  expression  for  oxidizer-fuel 
gas  reaction  distance 
Oxidizer  monopropellant  gas  reaction 
distance 

Oxidizer  gas  column  radius 
Oxidizer-Fuel  gas  reaction  distance 
Average  surface  mass  flux  (based  on 
planar  area)  at  a given  increment 
mass  flux  of  fuel 


“ox  ,p 


"‘ox , s 


MU  fuel 


MW„ 


oxidizer  mass  flux,  based  on  planar  sur- 
face projection 

oxidizer  mass  flux,  based  on  actual  total 
curved  surface  area 

molecular  weight  of  fuel  gases  leaving 
propel lanl  surface 

Molecular  weight  of  oxidizer  gases  leav- 


ing propellant  surface 
Global  gas-phase  reaction  order 
Pressure 

Total  heat  flow  to  surface  from  gas-phase 
reactions  (energy/time) 

^diffusion  flame  Beat  flux  to  surface  from  fuel- 
oxidizer  gas  flame  (energy/area/time) 
flame  Beat  flux  to  surface  from  com- 

bined monopropellant  and  fuel-oxidizer 
gas  flames  (energy/area/time) 

Heat  release  via  subsurface  reactions 
(energy /time) 

Heat  release  associated  with  HClO^(g)  + 
NH3(g)-  Equilibrium  Products 


P 

^gas 


^series 


‘I  sub 

Qap 


QeXO  heat  release  per  unit  mass  of  oxidizer 

consumed  in  surface/subsurface  reactions 
Qmelt.f  heat  of  melting  of  binder 

Qmelt.ox  heat  of  melting  of  oxidizer 

Qfuel.vap  heat  of  pyrolysis  of  fuel  binder 
Qsubl  heat  of  sublimation  of  oxidizer 
R(j  fraction  of  AP  reacted  per  unit  time  in 

subsurface  reactions 

Rgs  outer  radius  of  fuel  annular  gas  column 

rfuel  linear  regression  rate  of  fuel  surface 

rox  linear  regression  rate  of  oxidizer,  normal 

to  its  surface 

SMRBO  stoichiometric  ratio  (mass)  of  binder  to 
oxidizer 

T temperature 

Tf  flame  temperature 

T0  propellant  bulk  temperature 

Ts  surface  temperature 

TAU  time  for  fuel  regression  plane  to  move 

from  one  increment  to  the  next 
^crossflow  mainstream  crossflow  velocity 

Ecrossflow.x-  FH90  sin  0 

linear  gas  crossflow  rate  at  distance 
FH90  sin  0 from  the  surface 


Bgas,x-FH90  sin  0 

linear  gas  flow  rate  away  from  propellant 
at  distance  FII90  sin  0 from  the  surface 

U dimensionless  crossflow  velocity  at  x -- 

FH90  sin  0 

vgas  ,surf  gas  velocity  away  from  propellant  surface 

VLO  volumetric  fraction  of  oxidizer  in  pro- 

pellant 

WPO  weight  fraction  of  oxidizer  in  propellant 

X distance  from  surface  (sign  convention 

such  that  it  is  negative  below  the  surface) 

XDTOP  distance  of  fuel  surface  from  initial  peak 
of  oxidizer  particle  (See  Fig.  1). 

V*  dimensionless  value  of  FH90  sin  0 

Y2,o  mole  fraction  of  fuel  in  binder-pyrolysis 

products  in  absence  of  subsurface  react- 
ions. 

Yl,0  net  ' °le  fraction  of  oxidizer  in  oxidizer 

decomposition  products  in  absence  of  sub- 
surface reactions  (corrected  for  fuel  mole 
fraction  initially  in  the  oxidizer  stream, 
e.g.,  NHj,  which  negates  part  of  the  oxi- 
dizing value) 

-»  mass  fraction  of  oxidizer  which  reacts  at 

or  below  the  surface 

P mass  fraction  of  binder  which  reacts  at 

or  below  the  surface 

0 angle  between  resultant  velocity  vector 

and  planar  surface  (See  Fig.  2). 


a 


■"futtl 


\ gas  thermal  conductivity 

Xox  solid  oxidizer  thermal  conductivity 

P[  fuel  (binder)  density 

0OX  oxidizer  density 

“propellant  Propellant  density 

♦ (actual  molar  fuel/oxidizer  ratio) /(stoi- 

chiometric molar  fuel/oxidizer  ratio) 
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Figure  3.  Predicted  and  Experimental  Burning  Rates  for  Formulation  4686 
( 73/27  AP/HTPB,  5 Micron  Diameter  AP  ) Using  Initial 
Estimates  for  Subsurface  Reaction  Constants. 
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Figure  6.  Predicted  and  Experimental  Burning  Rates  for  Formulation  4685 
( 73/27  AP/HTPB.  5 Micron  Diameter  AP ) Using  Optimized 
Estimate  for  Subsurface  Reaction  Pre-exponential. 
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Figure  4.  Predicted  and  Experimental  Burning  Rates  for  Formulation  4525 
( 73/27  AP/HTPB,  20  Micron  Diameter  AP  ) Using  Initial 
Estimates  for  Subsurface  Reaction  Constants. 
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Figure  7.  Predicted  and  Experimental  Burning  Rates  for  Formulation  4525 
( 73/27  AP/HTPB  , 20  Micron  Diameter  AP  ) Using  Optimized 
Estimate  for  Subsurface  Reaction  Pre-exponential. 
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Figure  B.  Predicted  and  Experimental  Burning  Rates  for  Formulation  5061 
( 73/27  AP/HTPB , 200  Micron  Diameter  AP ) Using  Initial 
Estimates  for  tubiurfeea  Reaction  Constants. 
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Figure  8.  Predicted  and  Experimental  Burning  Rates  for  Formulation  5061 
{ 73/27  AP/HTPB,  200  Micron  Diameter  AP ) Using  Optimized 
Estimate  for  Subsurface  Reaction  Pre -exponential. 
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have  been  made  for  a series  of  three73/27  ammonium  perchlorate  (AP) /hydroxy- 
terminated  polybutadiene  (HTPB)  formulations,  with  oxidizer  particle  diameters 
of  S,  20,  and  200  microns,  and  compared  with  data  for  these  formulations.  With 
optimization  of  three  "free  constants"  appearing  in  the  model,  it  is  found  to 
give  excellent  agreement  with  no-crossflow  burning  rate  data  over  the  entire 
range  of  pressures  and  particle  sizes  studied.  In  all  cases,  however,  the 
predicted  sensitivity  of  burning  rate  to  crossflow  velocity  is  somewhat  less 
than  observed  experimentally. 
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